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Abstract 

Upper bounds for neutral as well as charged Higgs boson masses in a two-doublet 
model are obtained on the basis of tree unitarity conditions a la Lee, Quigg and Thacker. 
A wide variety of scattering processes are considered so extensively that our bounds are 
more restrictive than those obtained previously for neutral Higgs bosons and are also of 
a new kind for charged Higgs boson. It is argued that at least one of the Higgs bosons 
should be lighter than 580 GeV/c 2 . 
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1. Introduction 

Although the success of the standard SU(2)w x U(l)y gauge theory of the electroweak 
interactions is overwhelming, the Higgs boson sector in charge of spontaneous symmetry 
breaking has so far eluded experimental verification and is still a mystery. We all agree 
that the Higgs boson is one of our central concerns of the present experimental search and 
will be so even more in the future colliders, JLC, SSC and LHC. 

It has been known by now rather well that the mass of the Higgs boson, which is 
proportional to the Higgs quartic coupling, may be bounded from above, provided that 
the quartic coupling is not so large as to violate the validity of perturbative calculations 
[1,2]. In fact in the minimal standard model with a single Higgs doublet, Lee, Quigg and 
Thacker (LQT) [1] deduced the constraint from the perturbative unitarity which turned 
out to be m h < (8nV2/3GF) 1/2 = M L q T ~ ITeV/c 2 , where G F is the Fermi constant 
and mh is the mass of a neutral Higgs boson. 

Extension of their type of analyses was considered by several authors in the presence 
of more than one Higgs doublet. There are several motivations for increasing the number 
of Higgs bosons; supersymmetric extension of the standard model, a model of spontaneous 
CP-violation due to Higgs sector [3], the Peccei-Quinn mechanism [4] and so forth. There 
is also a surge of phenomeno logical interest in the two Higgs doublet model in recent 
literatures [5]. 

Casalbuoni et al. [6] have raised a question which has close bearings on LQT's, 
namely, at what energy strong interaction phenomena would start to show up whenever 
one or more Higgs masses are sufficiently large. They examined models with two doublets, 
a doublet plus a singlet and also a supersymmetric model where there exist three Higgs 
supermultiplets . 

Maalampi et al. [7] have recently studied the two-doublet model in the same vein 
as LQT. They derived an upper bound of the neutral Higgs boson mass by a numerical 
analysis, which gave them more or less the same bound as of LQT. It should be pointed 
out herewith that they did not consider a broad class of scattering processes to derive 
constraints on all of the charged and neutral Higgs boson masses. 

The purpose of the present paper is to reexamine the two-doublet model to see whether 
one can derive an upper bounds for neutral as well as charged Higgs boson masses in the 
method of LQT. We will answer to this question in the affirmative by taking into our 
considerations sufficiently large number of scattering processes. Overall allowed regions 
of three neutral and one charged Higgs boson masses are explored and their maximally 
possible values are presented (see Eqs. (45)-(48)). Moreover we will argue that at least 
one of the Higgs bosons ought to be lighter considerably than might have been expected 
from the LQT' work (see Eqs. (49) and (50)). 



2. Two Higgs Doublet Model 

Let us start by specifying the SU(2)w x U(l)y invariant Higgs potential for two 
Y = 1, SU(2) W doublets, $i and $ 2 - To avoid the flavor changing neutral current, we 
assume the discrete symmetry under $ 2 — ► — $2 [8]. The most general potential then 
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consists of five quartic couplings together with mass terms 
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V($ u $ 2 ) = XX-^I^I 2 + A ^l^| 4 ) + A 3 |$i| 2 |$ 2 | 2 
i=i 

+ A 4 (Re$J$ 2 ) 2 + A 5 (Inr$f$ 2 ) 2 . (f) 

The spontaneous symmetry breaking is triggered by two vacuum expectation values, v\ 
and i>2 of each doublet field and we write 

*i= (j^,.rt-L,v/>V ( 2 ) 



^-(Vi + /l< + IZi) 

In general, we can take i>i > 0, and v 2 > 0. 

The mass terms in (1) may be diagonalized by rotation 

hi \ _ / cos a — sin a \ / h \ 
hi ) \ sin a cos a y \H J ' 

wi \ _ / cos /3 — sin fi\ f w\ ( z i \ _ ( cos — sin (3\ f z 
w 2 )~\ sin/3 cos/3 J V G J ' V z 2 / V sin/3 cos/3 J \£ J ' 

In fact the mixing angles are determined by 



(3) 
(4) 



-Ai^ 2 + A 2 f| + aA^ 2 - A 2 ^)2 + (A 3 + A 4 ) 2 ^i rz\ 

tana = — — — — . {o) 

(A 3 + M)viv 2 

(— 7r/2 < a < 7r/2) and tan/3 = V2/V1 (0 < /3 < 7r/2). We note that w and 2; are 
massless Nambu-Goldstone bosons and are absorbed into the longitudinal components of 
the gauge bosons. The masses of the other fields, h, H, G, and £ are expressed in terms 
of Xi (i = 1, • • • , 5) and Vi (i = 1,2). Or conversely, the five quartic couplings are given 
by [6] 

. G JP ( 9 9 2 • 2 1 / r>\ 

Ai =— -j= m ft cos a + m H sin cm, (6) 



. Gf / 2 ■ 2 2 2 I — 

A 2 = ; 2 — I m h sin a + m H cos a], (7J 



\/2cos 2 /3 
V2sin 2 /3 

A3 = ^g 2 " (mj - mi) + 2V2 Gf „4, (8) 
A 4 = - 2V2G F m 2 G , (9) 
A 5 =2y/2G F (ml - m^j , (10) 



where we have set a/i> \ + v \ = (V2Gp) and have assumed j3 ^ 0, 7r/2 v 2 > 0). 



We will take it for granted hereafter that the quartic couplings Aj (i = 1, • • • 5) are 
lying in the region ensuring the positivity of all the masses squared. The conditions are 
summarized by 

4AiA 2 > (A3 + A4) 2 , (11) 
Aiu? + \ 2 vl > 0, A 5 > A 4 , > A 4 . (12) 
We will also assume rrih > ran without spoiling generality. 



3. Eigenvalues of S-Matrix 

The formulae (6) - (10) show that, if the quartic couplings \ (i = l,---,5) are 

— 1/2 

comparable to, say, unity, then the Higgs masses are also on the order of G F . The 
upper bounds of neutral as well as charged Higgs boson masses are derived by assuming 
that all these quartic couplings are within a perturbative region. As a criterion of the weak 
quartic couplings we adopt the perturbative unitarity a la LQT, thereby constraining the 
diagonal elements of the S-matrix. 

It has been known that a potential threat to unitarity comes from longitudinal gauge 
bosons and Higgs particles. Scatterings involving the longitiudinal gauge bosons are re- 
placed in the high energy limit by those of correponding Nambu-Goldstone bosons thanks 
to the equivalence theorem [1, 9, 10]. Our problem therefore boils down just to focusing 
the Higgs- Goldstone system described by (1). Furthermore, the dominant contribution in 
the high energy comes only from the quartic couplings, which we will consider henceforth 
exclusively. The S-matrix then becomes independent of energy, i.e., just a constant. 

We are interested in various two-body scatterings between fourteen neutral states, 
w + w~ , w + G~, G + w~ , G + G~, zz, zQ, hh, hH, HH, hz, Hz, Qh, and (H. The 
interactions in (1) are, however, extremely involved in terms of w, G, z, £ , h, and H, 
and so will be the S-matrix, too, not to mention the diagonalization thereof. It is now 
important to notice that the S-matrix evaluated in the mass eigenstate bases can always be 
transferred into the one in the original fields Wi, Zi, and hi (i = 1, 2) by making a unitary 
transformation. The S-matrix in the original field bases is to be calculated by using the 
quartic part of the interactions (1) expressed in terms of the original fields which looks 
much simpler of course. Since all we need to know are the eigenvalues of the S-matrix, it 
suffices to deal with the S-matrix in the original field bases. We are thus able to streamline 
our calculations just by using the original fields, Wi, Zi, and hi (i = 1, 2) from the outset. 

The use of the original fields are justified by the following considerations. In the high 
energy scattering, dominant contributions to the amplitudes come merely from the quartic 
coplings as was mentioned before, i.e., Feynman diagrams containing triple Higgs couplings 
are suppressed in energy on the dimensional account. The absence of Higgs propagators 
in our calculations amounts to ignoring admittedly the mass differences of various Higgs 
states originating from propagators. As in the usual particle mixing among degenerate 
states (such as massless neutrinos), we are allowed to choose freely the most suitable field 
variables. 

We evaluate scattering matrix on the basis of the states, \wfw~ >, \ziZj >, \hihj >, 
\hiZj > (i,j = 1,2) and/or combinations thereof. We will find following orthonormal 
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bases the most convenient: 



M >= 



Bi >= 



-^j=\2wfw i + ZiZi + hihi >, 



1 



2V2 



2wfw] 



1 



v/2' 



F± >= -\wiW 2 +W2 , w 1 ±ziz 2 ±hih 2 >, 

F x >= -k=\z\z 2 - h\h 2 >, \F 2 >= -^=\hiz 2 + h 2 z x > . 
V2 V2 



(13) 

(14) 

(15) 
(16) 
(17) 

(18) 
(19) 



The tree level high energy scatterings consist dominantly of S-wave, since they are 
described by the contact interactions with quartic couplings in the above bases. The transi- 
tion matrix of S-wave amplitudes turns out to be of a block-diagonal form diag (A, £?, C, D, E, f + , /_, fi, f 
with the following submatrices 



1 



6A1 



16tt V2A3 + i(A 4 + A 5 ) 



2A 3 + ±(A 4 + A 5 ) 
6A 2 



B 



16tt 



C = D = 



2Ai |(A 4 + A 5 ) 
KA4 + A5) 2A 2 



2Ai o(^4 — A5) 



16tt l^(A 4 - A 5 ) " 2A 2 



E 



1 /A 3 + ±(3A 5 -A 4 ) 



16tt 



-A, 



-A, 



A 3 + ±(3A 5 -A 4 ) 



(20) 

(21) 
(22) 
(23) 



The rows and colums of these 2x2 matrices are spanned by the states, (13)-(17), respec- 
tively. The diagonal elements of the amplitudes /±, fi and f 2 are those of (18) and (19) 
and are given by 



/+ 



^(As + ^-i^), /_ = ^(A 3+ iA 4 -iA 5 ), 
h = h = T^3+^ 4 + iA 5 ). 



(24) 



(25) 



We are lucky enough to obtain the fourteen eigenvalues of the transition matrix analytically. 
The eigenvalues of (20)-(23) are given respectively as follows; 



a± 



T ^{3(A 1 + A 2 ) ± y9(A 1 -A 2 ) 2 + {2A 3 + ^(A 4 + A 5 )} 2 }, 



(26) 



f>± =i^{(Ai + A 2 ) ± \/(Ai - A 2 )2 + i(A 4 + A 5 ) 2 }, (27) 



c± = <i± =7^{(Ai + ^) ± \/(Ai - A 2 ) 2 + i(A 4 - A 5 ) 2 |, (28) 

ei= lk (M - 1 2 ) " +5 2^' (29) 

e2= li (A3 -^ + I A5) ' (30) 

The reason for the block diagonal form of our transition matrix is explained in the 
following way. First of all let us recall the discrete symmetry under $ 2 — ► — $2- This 
indicates that the number of suffix "2" is conserved modulo two throughout scattering 
processes. In other words, the states \Ai >, \Bi >, |Cj >, and \Di > are not able to 
communicate with the other states \Ei >, \F± >, and \F^ >. 

Besides this discrete symmetry the quartic part of Higgs potential (1) is endowed with 
the following discrete symmetry 

4 (31) 
exp(^7r)^, (32) 

exp(^7ra 2 )$J. (33) 

The various states (13)-(19) are classified according to the properties under these discrete 
transformations. It is straightforward to see that \Ai >, \Bi >, \Ci >, \F± > and |Fi > 
are C-even states, while the others are C-odd. We also note that \Ai >, \Bi >, \Ei > and 
\F± > are Y^-even and the others are Y^-odd. On the other hand, not all of the states 
are the eigenstates of the G-transformation: there are only two G-even states (\Ai > and 
\F+ >) and two G-odd states (\Bi > and \F- >). These discrete symmetries are efficient 
enough to set up selection rules and to understand vanishing elements of the transition 
matrix between those states with different quantum numbers. 

In their analysis of the S- matrix, LQT argued that there was an 0(4) symmetry in the 
single-doublet Higgs potential in the high energy limit. They made a full use of this in close 
analogy with isospin symmetry, classifying states into 0(4) representations. Although our 
states (13)-(16) are reminiscent of their (9(4)-analysis, our Higgs potential (1) does not 
have such continuous symmetries. 

In passing our S-matrix was once studied partly by Casalubuoni et al. [6]. They 

— 1/2 

restricted themselves, however, to a special case that was much greater than G F 
and the mixing angle a was negligibly small. Note also that Maalampi et al. [7] have 
studied scattering processes with seven elastic channels w + w~ , G + G~ , zz, hh, hH, 
and HH (diagonal channel analysis). 



C : 




Y n : 




G : 





4. Analysis of the Tree Unitarity Constraints 
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Now that we have the eigenvalues of the transition matrix, we are in a position to 
impose the tree unitarity conditions . The unitarity is respected if the eigenvalues are 
lying within the so-called Argand circle on the complex plane. The tree level amplitudes, 
however, are necessarily real and therefore outside the circle. If higher order radiative 
corrections are included, the amplitudes will eventually settle down within the circle and 
unitarity will be restored. For such perturbative calculations to be meaningful, the tree 
level amplitudes should not be very far away from the circle. As a criterion for the per- 
turbative recovery of unitarity, LQT argued that the eigenvalues should not exceed unity, 

|a±|, |6±|, |c±|, |rf±|, N, |/±|, \fi\ < 1. (34) 

In some literatures [11], the RHS of (34) is replaced by 1/2 , the radius of the Argand circle. 
We, however, use more conservative conditions (34), LQT's, since it will be instructive to 
see how the inclusion of extra Higgs bosons will modify the LQT's result. 

For illustration, we take up the first one |a+| < 1. If we express Aj's in a+ in terms of 
nT>h, ma, ma, m^, a, j3, and G F , with the help of Eqs. (6)-(10), this inequality provides 
us with a relation that must be satisfied by various Higgs bosons masses together with the 
mixing angles. After a little manipulation we find \a + \ < 1 equivalent to the following two 
conditions 

9 ' x _M 2 _(V_y y> B 2, (35) 



9-5 sin 2 2a \ 3 
— sin 2 2/3 > X - Y cos 2a cos 2/3, (36) 



where our notations are 



X=-LG F (m 2 h + m 2 H ), (37) 

Y =l=G F (ml-m 2 H ), (38) 

Y = = |247rcos2acos2/3 - V2G F (m 2 c + 2m 2 G ) sin 2a sin 2/3 1 , (39) 

9 — 5 sin 2a [ J 



R = o 1 167r sin 2a cos 2/3 H — %=G F (m 2 r + 2m 2 r) sin 2/3 cos 2a 1 . 

9-5sin 2 2a\ v 7 ^ J 



(40) 



The second inequality (36) excludes the right half of the hyperbola (35) on the (X, Y) 
plane and the allowed region is the shaded one in Fig. 1, where Y > (mh > mu) and 
X -Y = V2G F mjj > are taken into account. 

Fig. 1 shows clearly that the allowed region of m# and rrih is bounded for fixed values 
of mixing algles, and ma- We can read off bounds for ran and from the shaded 
region 

X - Y = V2G F m 2 H < X P , (41) 
X + Y = V2G F m 2 h < 2X Q , (42) 
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where X P and Xq are the X-coordinates of the points P and Q, respectively in Fig. 1: 



X P = ^ - I y/(9 - 5 sin 2 2a) (r o 2 + rJ) , (43) 



Xq 5 sin 2 2a 



24tt - (9 - 5 sin 2 2a) Y 



- \J(V - 5 sin 2 2a){5.R 2 sin 2 2a + (8tt - 3r ) 2 } 



(44) 



We have similarly analyzed other eigenvalues, and found after all that (41) and (42) 
are the most stringent conditions upon mu and rrih- Since the points P and Q are both 
on the left half of the hyperbola, Xp < 87r/3 and Xq < 8n/3 should hold in general. We 
thus have upper bounds for neutral Higgs boson masses 



m H < 



m h < y-^r) = M LQT , (46) 

whatsoever the values of a, (3, and mg. 

The relations (41) and (42) contain more information than (45) and (46). To elucidate 
this, let us have a closer look at (41). In the (Gfiti'h, G f titq, Gi??n 2 )-space, The inequality 
(41) tells us that the inner region surrounded by the solid curves in Fig. 2 are allowed 
exclusively. This means that the charged Higgs boson masses vtiq and are also bounded. 
In fact Fig. 2 shows apparently 



» G <(^) 1/2 = 4 A w (47) 

mc<(^) =VTiM LQT . (48) 

The most interesting bound is derived by a geometrical inspection of Fig. 2. We obtain 
the bound on the lightest Higgs boson mass among m#, mc and m^, 
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M(the lightest Higgs boson mass) < \ — M L q T 




LQTi (49) 



where x = vijv\ = tan/?. This is our main result. The upper limit in (49) is realized at 
point K on the surface in Fig. 2 where the vector OK is parallel to (1, 1, 1). 



5. Summary 

In the present paper we have investigated the consequences of applying the tree uni- 
tarity conditions to the two Higgs doublet model. We have seen that, by considering a 
wide class of scattering processes, not only neutral but also charged Higgs boson masses 
are bounded from above. Our results are listed in Eqs. (45)-(48). 

The most important results of our calculations is that at least one of the Higgs bosons 
should be much lighter than has been anticipated from LQT's work [1]. The maximum 
value of the RHS of (49) is reached when x = 1 (sin 2 2/3 = 1). It is therefore concluded 
that the lightest Higgs boson mass has to satisfy 

M(the lightest Higgs boson mass) < -^-=M LQT = 580 GeV/c 2 , (50) 

v3 

that is, the bound is of the one derived by LQT for a single doublet Higgs case. It is 

very likely that inclusion of more Higgs doublets would give us more tight upper bound. 

A comment to be made herewith is that, if we would use |a+| < 1/2 as a criterion 
of perturbative recovery of unitarity [11], all the bounds of Higgs masses that we have 
obtained would be reduced by a factor l/y/2. In particular, the upper bound of the 
lightest Higgs boson mass (50) would be replaced by M lqt /VQ = 410 GeV/c 2 . 

Our final remark is that, so far, we have assumed v\, vi ^ 0. If one of these vacuum 
expectation values e.g. V2 happens to vanish, then the neutral Higgs boson hi becomes 
massless at the tree level. This is a pseudo- Nambu-Goldstone boson because there does 
not exist any particular symmetry prohibitting hi from acquiring a mass on the loop level. 
Then it is expected that V2 becomes non-zero by the higher loop corrections as far as $2 
is coupled with $1. 
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FIGURE CAPTIONS 

Fig. 1 The allowed region of Higgs boson masses given by (35) and (36). Our variables are 
X = G F {m\ + m 2 H )/V2 and Y = G F (m 2 h - m 2 H )/V2. 

Fig. 2 Schematic view of the allowed region of Higgs boson masses run, and ma given 
by Eq. (41) for given values of mixing angles. The inner region surrounded by the 
solid curves is allowed. The vector OK is in the direction of (1, 1, 1). 
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